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ON ■ The optical manifold method to compute the one-loop effective action in a 

static space-time is extended from the massless scalar field to the Maxwell field 
in any Feynman-like covariant gauge. The method is applied to the case of the 
Rindler space obtaining the same results as the point-splitting procedure. The 
result is free from Rabat's surface terms which instead affect the ("-function 
or heat-kernel approaches working directly in the static manifold containing 
conical singularities. The relation between the optical method and the direct 
(-function approach on the Euclidean Rindler manifold is discussed both in 
the scalar and the photon cases. Problems with the thermodynamic self- 
consistency of the results obtained from the stress tensor in the case of the 
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INTRODUCTION 



In a recent paper, we have computed the one-loop thermal partition function of pho- 
tons in the Rindler wedge employing a local (^-function method directly in the Euclidean 
Rindler space. Although this approach produces thermodynamical quantities with the cor- 
rect high temperature behaviour requested by the statistical mechanics, the low temperature 
behaviour seems to remain different from that obtained with others methods. This can be 
seen by means of a direct comparison between the free energy following from the above 
cited approach and the same quantity obtained by the point-splitting renormalization pro- 
cedure for the stress tensor f||-§J. In particular, one sees that the direct (^-function approach 
gives, for the coefficient of the term proportional to T 2 , a result which is one-third of the 
point-splitting result. This discrepancy can be traced back to an identical discrepancy in 
the coefficients of the free energy of a minimally coupled massless scalar field propagating 
in the Rindler wedge 0-0]. 

It is important to remark that this problem does not arise from the particular method 
used in and [|IJ to compute the determinant of the small fluctuations operator which 
appears in the one-loop free energy. In fact, the same discrepancy has also been found in || 
using a completely different method to compute the determinant. Therefore, it seems to be 
intrinsic of the computations made directly in the Euclidean Rindler space. 

In the photon and graviton further drawback of the approach in [||] is the need 

of a more complicated regularization procedure due to the presence of gauge depending 
"surface" terms |§. Anyway, the results of improve previous results obtained using 
global heat-kernel approaches flQ,9j in the Rindler space, which is not able to reproduce the 
Planckian high temperature behaviour. 

There is another method which can be used to compute these one-loop quantities, and 



is the optical one |TTHl8|. In this approach, instead of computing the partition function 
directly in the static metric, one performs a conformal transformation in such a way that the 
resulting manifold has an ultrastatic metric. Then, one can compute the relevant quantities 
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in this "optical manifold" using heat-kernel, ^-function, or any other method and taking 
into account how the various quantities transform under conformal transformations. This 
method is particularly favorable in the Euclidean Rindler case, since this manifold has 
a conical singularity which can be quite tricky to deal with, whereas the related optical 
manifold has no singularity. However, there is more in this method than the mathematical 
content. In fact, it has been shown |l^JT^Jl7|JT8[1 that the canonical partition function of 



a quantum field in a curved background with a static metric is not directly related to the 
Euclidean path integral with periodic imaginary time in the static manifold, but rather it 



is equal to Euclidean path integral in the related optical manifold. In particular, in JT4 
it is shown that the statistical counting of states leads naturally to a formulation in the 
optical manifold. We can also notice that, as far as we know, the equivalence of the direct 
periodic imaginary time path integral formalism to the canonical formalism for computing 
finite temperature effects has been proved in ultrastatic manifolds only |19 



Therefore, the computation of the thermal partition function from Euclidean path inte- 
gral in the static manifold requires the knowledge of the Jacobian of the conformal transfor- 
mation. On regular manifolds this causes no trouble, since it is easily shown (see, e.g., PU| ) 
that the Jacobian affects only the temperature-independent part of the free energy. Instead, 
as we will see, when in the static metric there is a conical singularity the temperature de- 
pendence of the Jacobian could be less trivial and affect the temperature-dependent part of 
the free energy. 

Another important point is that the optical method produces thermodynamic quantities 
which agree with those obtained form the point-splitting procedure. This happens in the 
case of a massless scalar field conformally coupled in the Rindler wedge at least, but also, 
as we shall see, in the case of the photon field. 

In the first part of this paper we shall review the computation of the thermodynamical 
quantities of a massless scalar field in the Rindler wedge, comparing the point-splitting, 
local ^-function and optical results. In particular, we note that while the point-splitting 
approach can be applied for any coupling of the scalar field with the gravity, the optical 



approach is feasible only for the conformally coupled case, where it gives the same result 
as the point-splitting method. Moreover, the dependence on the coupling parameter, which 
disappears in the integrated physical quantities when the background is a regular manifold, 
in this case affects these quantities because of the presence of the conical singularity. On 
the other hand, the computations made directly in the static Euclidean Rindler manifold 
using the local ^-function technique is still limited to the minimally coupled case, and even 
in this case the result is different from the point-splitting one. 

In the second part of this work, we shall extend the optical manifold approach to the 
Maxwell field case. We shall show that there are two possible ways to do this, which are 
equivalent in the scalar case, but in the photon case could produce a different result. The 
difference of the two approaches is essentially in the definition of the gauge-fixing and ghost 
parts of the Lagrangian. In particular, we shall show that result expected counting the 
polarization states of the field can be obtained by defining the theory directly in the optical 
manifold. The other possibility is to define the partition function in the Euclidean Rindler 
wedge and only then perform the conformal transformation. We are not able to work out 
thoroughly this latter approach, because of a mathematical complication in the ghost action. 
Nevertheless, we argue that it gives a different coefficient of the term proportional to T 2 . 

In the final part, in addition to discussing the obtained results, we also show that the 
usual relations between thermodynamical quantities, which involve derivatives with respect 
to the temperature, lead to inconsistencies when applied to Rindler thermal states with 
temperature different from the Unruh-Hawking one. 

I. THE GENERAL PROBLEM IN THE CASE OF A MASSLESS SCALAR FIELD 

As we have said in the Introduction, the problem of discrepancy of the coefficient of T 2 
in the free energy appears already in the case of a massless scalar field propagating in the 
Rindler wedge. In fact, this problem seems to be independent of the field spin. Hence, we 
start discussing just this case, taking into consideration the point-splitting, local C-function, 
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and optical manifold approaches. 



A. Point-splitting approach 



We start considering the results produced by the point-splitting method. The point- 
splitting renormalized stress tensor reads (see, for example 0, continuing into the Rindler 
space the results obtained for the cosmic string) 

4 



1 



14407r 2 r 4 
+20(6£ - 



jj -1 ) diag(-3, 1,1,1) 
llll^V-lW, -1,1,1) 



(1) 



By integrating —^/gT®, we get a total energy which we shall compare with those following 
from the other methods: 
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Above, L y and L z are the (infinite) lengths of the transverse dimensions, and so A± = L y L z 
is the (infinite) area of the horizon. The parameter £ fixes the coupling of the scalar field with 
the gravitation. In the (Lorentzian) Rindler space the scalar curvature R is zero everywhere 
and the parameter £ remains as a relic of the fact that T^ u is obtained by varying the metric 
g^u in the field Lagrangian pifl . [] Employing the general expression of T^ u {^) P,|21"|j in terms 



l \t is worthwhile noticing that one has to consider the theory within the curved space time in 
order to discuss on the physics in the flat space-time. Anyhow, the extension of the theory to a 
curved space-time is not unique and this involves some subtleties regarding also the regularization 
procedure. The choice between different regularization procedures should be made on the basis of 
what is the physics that one is trying to describe. Obviously, the general hope is that, at the end 
of the complete renormalization procedure involving matter fields and gravity, all these different 
regularization approaches give rise to equivalent physical results. See [[17]] and [22| for a discussion 
on these topics. 



of the Hadamard function, one finds that, in the case R = 0, the global conserved quantities 
as total energy should not depend on the value of £. This is because the contributions 
to those quantities due to £ are discarded into boundary surface integrals which generally 
vanish. However, this is not the case dealing with the Rindler wedge because such integrals 
diverge therein.^ The only possibility to get a result not depending on £ consists in taking 
(3 = 2ti producing a trivial result. The considered ambiguity does not seem to arise from a 
similar ambiguity in defining the thermal quantum state. In fact, the thermal Wightman 
functions employed in calculating the renormalized stress tensor do not depend on £. Finally, 
it is worthwhile noting that the £ ambiguity affects the (3~ 2 term in the thermodynamical 
quantities and hence their low temperature behaviour. 



Notice that Kay and Studer p3| found an ambiguity in defining the scalar Wightman 
functions around a cosmic string, a background which has the same Wick-rotated metric as 
the thermal Euclidean Rindler manifold. However, this ambiguity is related to the time- 
independent modes and so, e. g., employing the ^-function procedure, one simply finds that 
this ambiguity cannot produce f3~ 2 terms in the Rindler free energy. Hence, it should not 
be related with the ^-ambiguity. 

B. Direct conical approach 

One can formally define the partition function at the temperature 1/(3 by an Euclidean 
path integral 

e ~^ = Zp = Jv<j>exp(-S[<l>]) , (3) 
where the (Euclidean) action is that of a massless scalar field coupled with the gravitation, 

S[<f>] = -\j d A x^ [N^V - £R] 0. (4) 



2 Similar problems appear working in subregions of the Minkowski space in the presence of bound- 
ary conditions [21]. 
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The background is the Euclidean Rindler manifold Cp x R 2 with an imaginary time period 
(3. The Euclidean Rindler metric reads 



ds z = r 2 d6 2 + dr z + dy z + dz< 



(5) 



where 9 G [0, f3], r G R + , x = (y, z) G R 2 . Notice the well-known conical singularity at r = 
when f3 7^ 2n. 

In the case £ = 0, the previous partition function can be explicitly computed by a 
local ^-function approach recently introduced by Zerbini et al. f| obtaining a Minkowski 
renormalized free energy F| ub = Fp — Up=2n and a renormalized internal energy^ ?7| ub 
dp(3Fp - {dp(3Fp)\p =2lT which read 
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where A± is the (infinite) event horizon area and e a short- distance cutoff representing the 



minimal distance from the horizon 241. We note that the above result is different from 



that of the point-splitting with £ = and the difference is in the coefficient of the term 
proportional to f3~ 2 . However, as we said in the Introduction, the partition function ([|) 
is related to the canonical one by the Jacobian of a conformal transformation: this could 
explain the different result. We shall come back on this point after discussing the optical 
approach. 

It is worthwhile noticing that the Lorentz section of the Rindler space is flat and hence, as 
far as the real time theory is concerned, we find a complete independence on the parameter 
£. However, in calculating the partition function, one has to deal with the Euclidean section 



3 As it is well known, the {(3 = 27r)-thermal Rindler state locally coincides with the Minkowski 
vacuum and, in renormalizing, we suppose that this state does not carries energy density. Notice 
that such a Minkowski subtraction procedure does not affect the entropy computed through Fp. 
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of the Rindler manifold and, considering it as a integral kernel, the curvature R takes Dirac's 
S behaviour at r = [^,^], thus the value of the parameter £ could be important. The 
previous results have been carried out in the case £ = in the sense that the eigenfunctions 
employed in computing the ( functions properly satisfy the eigenvalue equation with no R 
term. 

In the case £ 7^ the problems are due to the fact that the equation for the eigenfunctions 
contains a Dirac 5, and so it is not mathematically clear how to treat it. In the case of a 
cosmic string, the Dirac 5 represents a limit case, maybe unphysical, of the problem in 
which the string has a finite thickness, which is mathematically well defined since no Dirac 
S appears. In the case of the Rindler space there is no such way out, and the only way to 
avoid the problem is to consider the case £ = 0. 

C. Optical approach 



Let us now consider the optical approach [fTI-IB|. As remarked in the Introduction, 
this approach is not just a mathematical method to compute the functional integral in Eq. 

but has an important physical content. In fact, as previously remarked, it fulfills the 
requirements of a formulation in the optical manifold following from statistical counting of 
states. 

Let us consider a static metric ds 2 = g tiv dx^dx v and perform a conformal transformation 
of the metric (maybe singular if Q(x) = 0), g^ v — > g' = Q 2 (x)g lJiU , so that 

ds 2 -> ds' 2 = Q 2 ds 2 . (7) 

Choosing f2~ 2 = g 00 , ds' 2 becomes the related ultrastatic optical metric. In the case of 
the Euclidean Rindler space, this conformal factor becomes singular just on the conical 
singularities, which are pushed away to the infinity^ and the optical manifold is free from 



l The points at r = of the optical manifold S* 1 x are infinitely far from the points of the 
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singularities. Under such a transformation, the massless scalar field transforms into 0' = 
Vt~ l (f) and the Euclidean action with coupling factor £ transforms into the following more 



(8) 



complicated action Jzl 

s'W\ - 

where £b = (D — 2)/A(D — 1) is the conformal invariant factor. If we consider a conformally 
coupled field, £ = we see that also the transformed action is that of a conformally 
coupled field in the optical manifold S 1 x H 3 . In the other cases, we have to keep a term 
proportional to R which has a Dirac 5 behaviour at r = and thus we have an ill-defined 
operator. (] 

When we compute the one-loop partition function, we may formally write [27, 2q| 

e-W = Z P = J Vet)' J[g, g', (3] exp = J[g, g', (3]Z' p = J[g, g', (3]e~^ . (9) 



We remark that is Z' g which is equivalent to the canonical partition function Jl3|,|T^,[T^,|T8 



The functional Jacobian J[g, g', f3\ does not depend on <$' and thus it can be carried out from 
the integral as we have done above. When the involved manifolds are regular, it is possible 
to prove that such a Jacobian is in the form 



J[g,g',(3]=exp(-(3E ) 



(10) 



where E does not depend on (3. This is substantially due to the staticity of the involved 
metrics [16|-18 and the factor {3 in the exponent is due to an integration over the whole 



manifold with r > taking the distance as the affine parameter along geodesies. Strictly speaking, 
the former points do not belong to the manifold at all. 

5 One possible way to get rid of this term is to define the action in the Lorentzian manifold, 
where R = 0, perform the conformal transformation to the optical manifold, and only then use 
the transformed action to write to partition function with the periodic imaginary time formalism 
|15| . This procedure gives a result independent on the parameter £ by nature: the coupling in the 
optical manifold is always conformal. However, in our opinion this procedure seems too ad hoc. 
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Euclidean manifold. If this holds in the presence of a conical singularity as well, one expects 
that F@ and differ only for the value of the renormalized zero-temperature energy. When 
the coupling in the Euclidean Rindler manifold is conformal, the direct computation of F'p 
can be performed employing the (^-function approach ]I6| (see also the Appendix of this 
paper). We report here the well-known final result only: 

F' = — 
& 2880vr 2 e 2 ^ (3 

Using JJ'a = dp{(3F'p) to compute the internal energy and performing the Minkowski renor- 
malization L^ sub = Up — t/g =27r i n order to get a vanishing internal energy at f3 = 2n, we 
find just 

Therefore, we have got a result equal to the point-splitting one by renormalizing (with 
respect to the Minkowski vacuum) the internal energy obtained on the optical manifold and 
without taking into account the Jacobian, whether it has the form ( PD or not. 




D. Comparison of the results 

In the previous subsection we have seen that the optical method, when applicable, gives 
the same result as the point splitting. On the other hand, we see that Up'~^' does not coincide 
with the corresponding internal energy @ found by the ^-function approach at the value of 
coupling parameter one expects, £ = 0, but rather at £ = 1/9. Note that the discrepancy 
is in the term proportional to /3~ 2 , while the difference in the /^-independent term is not 
meaningful, because such terms are fixed by the subtraction procedure: they do coincide 
when the remaining terms are equal. Note also that we cannot compare directly the optical 
and the local ^-function approaches, since they are not applicable for the same value of £. 

In order to identify the source of the above discrepancy, we remind the reader that the 
thermodynamical quantities in Eq. (||) have been obtained from the Euclidean path integral 
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in the static manifold, which differs from the canonical partition function for the Jacobian of 
the conformal transformation, see Eq. (|9|). As we have said above, on regular manifolds the 
logarithm of this Jacobian is simply proportional to /3, thus giving a contribution only to the 
temperature-independent part of the free energy However, the case of the Euclidean Rindler 
space could be more complicated, due to the presence of a conical singularity at r = 0, which 
could yield a nonlinear dependence on (3. In fact, such a singularity can be represented as an 
opportune Dirac 5 function with a coefficient containing a factor (2ir — j3) p5j^6| , and so (3 
enters not only as integration interval, but also in the integrand. Of course, only an explicit 
calculation of the Jacobian can give an ultimate answer. In two dimensions, the Jacobian 
J[g,g',/3] is the exponential of the well-known Liouville action [pTfl , and an easy calculation 



shows that the logarithm of the Jacobian is indeed proportional to f3 ||18|| , regardless of 
the conical singularity. Unfortunately, in four dimensions the form of the Jacobian is far 
more complicated (see |18j and references therein) and involves also products of curvature 
tensors which are ill defined. Therefore, it is not clear whether the discrepancy in the term 
proportional to (3~ 2 might be assigned to the Jacobian. 

Summarizing, we have seen that the optical method has been applied to the conformally 
coupled case only, and in this case it gives the same result as the point-splitting method. 
With regard of the direct computation in the Euclidean Rindler wedge, it has to be con- 
sidered as incomplete, because of our ignorance of the Jacobian and of the nonminimally 
coupled case. We stress that, in the case of a regular manifold, these two approaches should 
be equivalent. 



II. OPTICAL APPROACHES IN THE CASE OF PHOTONS 

In J]]] the partition function of photons gas in a Rindler wedge has been computed 
generalizing the procedure in ||. The found Minkowski renormalized free energy amounts 
to 2F| ub + (2 — lna).F / ! urface , where F^ nh is the scalar free energy previously discussed, Eq. 
([]), and the _p| urface j s a "surface" term which arises integrating a total derivative and has 
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the form A±[(2n/(3) 2 — l]/(247r 2 e 2 ) (see |J and JT| for more comments), finally a is the 
gauge-fixing parameter. Notice that also this anomalous gauge-dependent term involves a 
j3~ 2 dependence. We suggested droping this latter gauge-dependent term as the simplest 
procedure to remove the unphysical gauge dependence. Anyway, we stressed that other 
procedures could also be possible. The obtained result agrees with the statistical mechanics 
request at high temperatures, but, as in the scalar case, the low temperature behaviour 
is different from that obtained with the usual point-splitting procedure. Therefore, let us 
consider the point-splitting results @||: the renormalized stress tensor takes a simple form 

1 
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diag(-3, 1,1,1). (13) 



The (Minkowski renormalized) internal energy corresponding to the previous photon stress 
tensor reads 

2tt\ 4 /2vr\ 2 



jjphot. p.-s. ° ^ 



p ~~ 14407T 2 e 2 

As far as the energy density is concerned, we have the following very simple relation: 



(14) 



(T o ) pho, P .-, = 2(r o ) p.-, (e = 0); (15) 

where on the right the stress tensor is that of a massless scalar field. It is worth while 
noticing that £ = takes place on the right-hand side instead of £ = 1/6. Hence, the 
energy density of the electromagnetic field does not amount to twice that of a conformally 
coupled scalar field, as one could naively expect considering that the electromagnetic field is 
conformal invariant in four dimensions. As far as the internal energy is concerned, we find 
the same unforeseen relation. However, as previously discussed, the integrated quantities 
should not have to depend on £ in more "regular" theories, restoring the naively expected 
relation between the considered quantities. 

Discussing the scalar case we have stressed the importance of the optical method: there- 
fore, now we go to investigate whether it is possible or not to get such an energy employing 
the optical-manifold method. There are two possible ways to implement this method. The 
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simplest one consists of defining the partition function directly as a functional integral on 
the optical manifold. However, there is another more complicated possibility: it consists of 
starting with a functional integral in the initial static manifold, performing the conformal 
transformation and finally dropping the functional Jacobian. This is, in fact, the simplest 
generalization of the results obtained in the scalar case. Both methods produce the same 
final functional integral in the simpler conformally coupled scalar case, but in the case of 
the Maxwell field the two procedures do not seem to be equivalent, as we shall see, due to 
the presence of gauge-fixing and ghost terms. 



A. Optical approaches in the case of general static manifolds 

Let us start reviewing the formalism we use dealing with the photon field. The complete 
action for the electromagnetic field in any covariant gauge and on a general Euclidean 
manifold, endowed with a metric ds 2 = g tiv dx ii dx v ', which we shall suppose static and where 
do is the global (Euclidean) timelike Killing vector with closed orbits of period (3. Using the 
in Hodge de Rham formalism we have 



S c 



d 4 x 



l(F,F) + ^-(A,d5A) 



ghost 



a 



d x 



(A,AA) 



1 



1 



a 



(A, d5A) 



+ S t 



ghost 



a 



(16) 



In order to maintain the gauge invariance of the theory, it is important to keep the de- 
pendence on the gauge-fixing parameter of the ghost action, as one obtains by varying the 
gauge-fixing condition -j^5A = |[29|| : 



S, 



ghost {Oi j 



—= [ d 4 x^fq cAc, 
JaJ v 



(17) 



where A is the Hodge-de Rham Laplacian for 0-forms and c, c are anticommuting scalar 
fields. Usually, the dependence on the gauge-fixing parameter is absorbed rescaling the 
ghost fields, but in the presence of a scale anomaly this rescaling gives rise to a nontrivial 
contribution, which is essential to maintain the gauge invariance of the theory. This is just 
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the case here: in fact, the contribution of the action fll7|) to the one- loop effective action is 
proportional to that of a minimally coupled scalar field, which has a scale anomaly in four 
dimensions. 

Some comments on the formalism in Eq. flOp are in order. F = d^A v — d v A^ = 
Vfj,A u — V U A^ is the 2-form representing the photon strength field, V M being the covariant 
derivative; the brackets stand for the p-forms Hodge local product: 

(G, H) = g a *h = Vag^-g^G^H^ 

For future reference we also define the internal product 

G-H = g^...g^G, 1 ...,H Vl .. Mp . 

We remind the reader that 5 = (— l) Ar ( p+1 ) +1 * d* is the formal adjoint of the operator d with 
respect to the scalar product of p-forms induced by the integration of the previous Hodge 
local product; finally, A = d5 + 5d is the Hodge-de Rham Laplacian of the p-forms. In 
order to perform calculations through the usual covariant derivative formalism the following 
relations for 0-forms and 1-forms are quite useful: 



A0 = -V^V^ 



5 A = -V M ^ , 



(AA)v = -V.VM, + R V A V . 



The second line of Eq. ([lq) represents the complete photon action now expressed in 
terms of the vector field A^ and the ghost fields only and it is the one usually employed in 
order to compute the partition function of the photon field by means of a functional integral. 
The partition function of photon at the temperature T = 1/(3 is then formally expressed by 



l4 x 



Z p = JvA exp--/ d 

x J VcDc exp -Sghost (a) • 



(A.AA) - (1-^) {A. (ISA) 



In order to compute this partition function, we want to pass to the related optical 
manifold, and so we consider a conformal transformation, Eq. (0), with Q 2 = goo- Notice 
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that, since we work in four dimensions, the p-forms A and F have a vanishing mass dimension 
and thus they must be conformally invariant, namely A = A' and F — F'. Furthermore the 
following identity arises: 

(F,F)'=(F,F). (18) 



B. First general approach 

As we said above, the way to proceed is twofold. As a first way, we can suppose to have 
performed the conformal transformation before we start with the field theory. This means 
that we define the partition function of photons in the static manifold as a path integral 
directly in the optical manifold. In such a case the expression of the partition function is 
defined by 

1 



Zf = J VA exp - X - J d 4 x (A, A' A)' - (1 - (A, dS'A) 1 
x J Vc'Vc exp-5^it(a) , 



(19) 



where 



y(i) 

ghost 



(a) 



1= ( d A x^'c'A'c' , 
s/a J v 



and where the primed metric and variables appearing in the previous functional integral are 
the optical ones. In other words, for the one-loop Euclidean effective action — In we 
have 



In 4" 



In det \i~ 



A' 



1 



a 



dS' 



+ ln 4!lhost(«) • 



(20) 



Here \x is an arbitrary renormalization scale necessary on a dimensional ground in the above 
formula and denoting the presence of a scale anomaly if it does not disappear from the final 
formulae. 

For future reference we note that the effective action of the ghosts, except for the a 
dependent factor, amounts trivially to minus twice the Euclidean effective action of an 
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uncharged massless scalar field with the Euclidean action minimally coupled with the gravi- 
tation. Therefore its contribution to the one-loop effective action can be written immediately 
from the ( function of a minimally coupled scalar field, C in,c,s, (s; x), in the same background, 
taking the a-dependence into account: 

d 



m 4!fko S t(") 



J d 4 x^J~g' 



ds 



C mx - s -(*; x) \ s=0 + C m - C - S is; x)\ s=0 In v^/i 2 



(21) 



C. Second general approach 

As a second way, we can suppose to define the partition function directly in the static 
manifold, adding also the gauge-fixing term and the ghost Lagrangian to the pure electro- 
magnetic action, and only after perform the conformal transformation to the optical metric. 
In this way we have to find how all the pieces in the path integral transform under the con- 
formal transformation. In particular, the operator A — (1 — a _1 )<i5 transforms into another 
operator A a , which we are going to write shortly. As regards the functional Jacobian which 
arises from the functional measure, a direct generalization of the discussion made in the 
scalar case tells us that it has to be ignored if we are interested in computing the thermal 
partition function. However, we would have to take it into account if we were computing, 
for example, the zero-temperature effective action in a cosmic string background. 

Hence, employing this second procedure, we shall assume the photon partition function 
to be defined by 



zf = J DA' exp- d 4 x(A',A a A')'\ J Vc'Vd exp -S' g 



? } (a) 

ghost I" I 



(2) 

In other words, for the Euclidean effective action — InZg we have 

InZf = -i lndet(/i- 2 A Q ) + \nZ^ host (a) . (22) 
The form of ^ghostW ^ s that of Eq. (|T7|) after a conformal transformation: 

5£ st (a) = ~ j d'x^'li [A' + 1{X - {l~ 2 R)\ c> , (23) 
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where d = Qc, c' = Qc. For future reference, we note that this effective action of the ghosts 
amounts trivially to minus twice the Euclidean effective action of an uncharged massless 
scalar field (p with the Euclidean action (A' = — V^V'^) endowed by an a— depending 
overall factor 



a 



A' + -(R' -n~ 2 R) 
6 



<P 



(24) 



When the static manifold is flat, R — 0, the contribution of the ghosts to the effective action 
can be written in terms of the ( function of a conformally coupled scalar field: 



ln 42ho S t(") = - / d 4 x^g' 



j- s C c - s is; x)\ s=0 + C c - S is; x)\ s=0 In v^u 2 



(25) 



Now, let us find the explicit form of the operator A a . The following identity holds: 

5A=±(5'A- V -A), 
provided the 1-form rj be defined as 



(26) 



7] = d(ln f2) = dfj, In 



(27) 



Taking into account that S = $ and employing Eq.s flTBD , (|18D and (|27|) we get the identity 



S 



"'" " 2 ' <Vx 

d^x 



(AAA)- 1 



a 



(A, d5A) 



Sghost(aO 



(A, A' A)' - (l--) (A,dd'Ay + -(A, VV -Ay 
V ay a 



— (A, ( V 5' + d V -)A)' 
a 



+ 5'gho ) st(«) 



(28) 



Looking at the first line of Eq. fl28|) we find the explicit form of the operator A a 

A„ = A' - (\ - dS' + -rjq ■ --(fjS' + drj-) . 
a J a a 



(29) 



Notice that the use of such an operator is equivalent to employing an unusual gauge-fixing 
term in the initial photon Lagrangian which reads 

1 



a 



■(A, (d-r})(5'-r}-)A). 



(30) 
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III. THE CASE OF THE RINDLER SPACE 



Let us check the physical results arising from Eq.s (^0) and Q2"2| ) in the case of the Rindler 
space. Setting Q 2 = r 2 in Eq. (0), the related ultrastatic optical metric reads 

ds' 2 = dT 2 + r~ 2 (dr 2 + dy 2 + dz 2 ). (31) 

Obviously, this is the natural metric of S 1 x if 3 which does not contain conical singularities. 
We remind one that R'Jf = —2 diag(0, 1, 1, 1) and R' = —6. As for the 1-form 7] necessary 
to define the operator A a , we get 

V, = ■ (32) 



We want to employ a local (^-function [30,31,20] regularization technique and hence we 



define the determinant of an (at least) symmetric operator L through 

1 1 f /— r i 

- -lndet(/i~ 2 L) = - / d A x^g' £(s = 0; x) + £(s = 0; x) ln/x 2 , (33) 

where the local ( function of the operator L is defined, as usual, by means of the analytic 
continuation in the variable s G C of the spectral representation of the complex power of 
the operator L: 

((s;x) = J2K s An(x)-A* n (x). (34) 



Above, A n (x) is a 1-form eigenfunction of a suitable self-adjoint extension of the operator 
L and A n is its eigenvalue. The index n stands for all the quantum numbers, discrete or 
continuous, needed to specify the spectrum. The set of these modes is supposed complete 
and (Dirac, Kroneker) 5 normalized. We will make also use of the following notation for the 
1-forms on S 1 x H 3 : 



A = (a\B) 



where a indicates a 1-form on S 1 and B a 1-form on H 3 . All the operations between forms 
which appear after "|" are referred to the manifold H 3 and its metrical structure only. Latin 
indices a, b,c,d, ... are referred to the coordinates r, y, z on H 3 only. 
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A suitable set of eigenfunctions of the operator A' — (1 — a" 1 )^' as well as A a as can be 
constructed using the following complete and normalized set of eigenfunction of the scalar 
Hodge de Rham Laplacian on S 1 x H 3 : 

„tkx „iu„T 

0( k <"^(r,r,x) = __^sinh(7ru;) rK ioJ (kr) , (35) 

where v n = ^zp, n G Z, u G R + , k = (k y ,k z ) G R 2 , k = |k| and all the previous eigen- 
functions have eigenvalue [y 2 + u 2 + 1). K iu) (x) is the usual MacDonald function with an 
imaginary index. The normalization reads 

J d A X^> 0(M,")*0(k',nV) = 5 nn> ^2 (k _ _ ^ 

In the following, we report some relations which are very useful in checking the results 
which we shall report shortly. It is convenient to define the 1-form £ = —d(l/r) = i]/2r 
on H 3 . On H 3 we have: 5r] = 4, Ar] = 0, drj = 0, rf£ = 0, A£ = -3£, V a £ 6 = -8%/r. 
Furthermore remind that, if / is a 0-form and u an 1-form: 

[A(/o;)] a = f[Au] a + u a Af - 2(V b f)V b u a . (36) 

Finally, on a 3-manifold the following relation holds 

A * (u A J) = *[(Aw) A «/] + A Aa;'] + i? * (u A «/) - *[(Riu) A u/] 

- * (u A fk/) - 2 * (V a cu A VV) , (37) 

where obviously [*(V a ^ A V a a/)] e := y/ge e b c (V dUJ b )'V d uj' c , u and a/ are 1-forms and the Ricci 
tensor acts on 1-forms trivially as (Ruj) a = R b a Ub- 

A. First optical approach 

Let us now consider the first optical approach, in which we define the path integral 
directly in the optical manifold, see Eq. (fK)|). Starting from the scalar eigenfunctions, one 
can obtain the following set of eigenfunctions of the operator A — (1 — a~ r )d5 on S 1 x H 3 . 
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\v\\Juj 2 + V 2 + 1 

1 



( d T (j) | d(j) 



\v\\Juj 2 + V 2 + 1 



(<9 T 0, d r (f>, d y (f), d z 4>) 



Vuj 2 + V 2 + 1 



Sir 



CJ 2 + 1 



^ = ^(0| *^0)) = ^(O,O,^,-^) 

AW = J-( | 5dm ) = -^(0, -0, d r d/, 

/ecu /ecu r r r 

The last three modes are transverse, 5A = 0, whereas the first one is a pure gauge mode. 
From a little Hodge algebra, the following normalization relations can be proved: 



d A x{A^ n * ] \ A (JV ' n '' k,) ) = 5 JJ '5 nn ' 5 2 (k - k') <f(w - J). 



(38) 



As far as the eigenvalues are concerned, we have: 

LU 2 + V 2 + 1 



A' - (1 - cT 1 )^' v4 (1) 



Of. 



A' - (1 - oT 1 )^'] A< 2 > = (cu 2 + z/ 2 + l)^ 2 ) , 

A' - (I - a~ l )d5'} A {J) = (v 2 + lu 2 )A {j) if J = 3,4 . 



Employing the definition in Eq. fl3~4), the above modes and the definitions given in the 
appendix, we have that (notice that <ft* and <ft take the same values of k, n, u) 

<j>*(x)<j>(x) ~ f j2 , J 2(l + u- 2 )<f ) *(x)<f ) (x) 



C(s;x) 



la 



+ 1) X) / ^krfcu 



[cj 2 + v 2 + !]•' 



+ ]C / d 2 kdoo- 



[uj 2 + u 



2 Is 



( a * + l)C m - c - s -( s; z ) + 2C c ' c - s -(s; x) + C extra (s; x) 



(39) 



where we have set 



( cxt ™(s; x) = 2 X / rf2k / 



gL> (jf(x)(j)(x) 
[u 2 + v 2 \ s 



2s-l 



" r(S -^#| 6,(2.-1), 



(40) 



7T 2 /? r(s) v 27r 

so that ( cxtra ( s = 0; x) = and ( /extra ( s = 0; x) = 1/3/3 2 . Notice that the second and third 
terms in Eq. ( |3"9"D arise from the transverse modes A^ and A^\ The first term in Eq. 
is due to the modes with J = 1, 2. 
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In calculating Eq. fl39|), we encountered Rabat's surface terms similar to those we en- 
countered in However, in the present case all these terms vanish automatically and no 
further regularization procedure needs. In fact, all these terms read as 



D r ^2jdkJ dur 2 K iuJ (kr)K iuJ (kr)f(uj,v, s) 



where D r is an opportune differential operator in r. Passing from the integration variable k 
to the integration variable rk, we see that the term after the operator does not depend on 
r, and so the differentiation produces a vanishing result. 

In order to write the complete local ( functions of the electromagnetic field we have to 
take account of the ghost contribution. We have already said that in this approach the £ 
function of the ghosts is just minus two times the ( function of a minimally coupled scalar 
field, but with a gauge- fixing dependent scale factor, see Eq. (|2T|): 



Using this relation, Eq. (^Uj) and reintroducing everywhere the renormalization scale /i, we 
can write the complete local ( function of the electromagnetic field as 



C m (s;x) = (a s + l)r c - s is^- 2 L^ )(x) + 2C c - s is^- 2 L^)(x) 

+C extra (s;/x- 2 )(x) - 2C m ' c - s -(s;/i- 2 «-^ 5 =o)(x). (41) 



It follows that the one-loop effective Lagrangian density is just 

£ cS (x) = l^\2C- c - s is-,x)+C xt ™(s;x) 

ids L J s=o 

= — + — (42) 

We remark the importance of keeping the a dependence of the action of the ghosts: it gives 
a contribution proportional to In a which cancels against the (lna)-dependent term coming 
from (a s + l)£ mx ' s '(s; x), restoring the gauge invariance of the theory. Note also how all the 
terms containing In fi 2 cancel giving the expected scale invariant theory. 
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Integrating this quantity over the manifold and introducing a cutoff at a distance e form 
r = in order to control the horizon divergence, we get the one-loop free energy: 

a r /~ \ 4 fn \ 21 

F^ = --Jd A x^'C eS (x) 



.7 + » 7 



(43) 



/?./ ~ ~v*~-w 14407r 2 e 2 

Renormalizing this result in such a way that the internal energy vanishes at /3 = 2tt, we find 
just the point-splitting result 



T7-(l)sub rrphot. p.- 

Up —Up 



(44) 



B. Second optical approach 

Let us then consider the second optical approach. We were able to perform the cal- 
culations in the case a = 1 only, hence a complete discussion on the gauge invariance (a 
invariance) is not possible. However, the found result contains some interest. As before, 
the eigenf unctions of the operator A a= \ are constructed from the scalar eigenf unctions, Eq. 

a® = ^(d T 1 + M#) = -L(c> T 0,|H-,o,o) 

\u\ 2 \u\ r 

A® = ±-Ad T <p\ - M#) = i-^-H-^^) 

\v\ 2 \v\ r 

k k r r 

A W = I( | * (£ A *d(^)) ) = i(o, o, a*-, d£) . 

k k r r 

The following normalization relations hold: 

J d i x{A^ n ^\ A( J '^'' n '' k ')) = 5 JJ '6 nn ' 5 2 (k - k') 5(lo - lo'). (45) 

As far as the eigenvalues are concerned, we have: 

A a=1 A^ = {u 2 + [(-l) J + \v\] 2 }A^ if J = 1,2, 
A a=l A^ = (z/ 2 + u 2 )A^ if J = 3, 4 . 
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Employing the definition in Eq. ([34]) and the found modes we have (notice that <fi* and 
take the same values of k, n, to) 

(x)<P(x) , ~ t r 20*(x)0(x) 



n=l 

oo 



+ + !)>]. j [uJ 2 + {u _ m 



[ . 40*(x)0(x) 



n=l 



[W 2 + Z/ 2 ] s 

For simplicity, we have omitted the terms corresponding to n — 0, which contribute only 
to the temperature-independent part of the free energy: this part will be changed during 
the renormalization process (subtraction of the Minkowski vacuum energy). We also stress 
that Rabat's surface terms involved during the calculations disappeared exactly as in the 
previous approach. The latter term in Eq. (|39| ) is due to the modes with J = 3, 4: this term 
is exactly twice the ( function of a conformally coupled Euclidean scalar field propagating 
in S 1 x H 3 . 

As far as the ghost contribution is concerned, it arises from the action (f23"D. Since the 
corresponding small fluctuations operator involves the curvature of the Euclidean Rindler 
manifold, which has a Dirac 5 singularity at r = 0, mathematically it is not well defined and 
is not clear how to deal with it. However, as a try we can suppose to consider R = and see 
the consequences.^ Under this hypothesis, the ghost contribution is just minus twice that of 
a conformally coupled scalar field (see Eq. (p5|) ) and so it cancels against the contribution 
of the modes J = 3, 4. 

After having added the ghost contribution, we can write the complete £ function of the 
electromagnetic field as 

<-(«;*) = E J <* J rf v + (, + 1) V + E J dk J + {v _ m . ■ < 47 > 

The partition function of the photons is obtained employing the previous function oppor- 
tunely continued in the variable s in Eq. (0). Dealing with it as in the previous case, we 
finally find the free energy 



6 See footnote number 5. 
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'2n 

J. 



-30 



2tt 



29 



(48) 



In deriving this result we have employed the Riemann zeta function ((z, q) and its relation 
with the Bernoulli polynomials ||32|| . This result has the same form as that obtained with 
the first approach, Eq. (fj5|), but the sign in front to the second term is opposite. The third 
term is fixed by the renormalization procedure. The problems arise with the f3~ 2 term once 
again. 

In this case it is easy to identify the origin of the discrepancy in our hypothesis of setting 
R = in the ghost action. Nevertheless, there is some evidence that the origin is not that. In 
particular, if we assume that the optical method gives the same results as the point-splitting 
one even when £ ^ 1/6, then we can suppose that it is right to substitute the optical result 
for the ghost contribution to the above free energy with the point-splitting one for £ = 0. 
As a result, we get 



_p(2)- 



sub 



1440vr 2 e 2 



2tt 



2tt 



— - 60 — +59 







ft 



(49) 



which is different from the previous one but still different from the first optical approach 
result. In particular, the free energy in Eq. ( [4"8"D (or (fj5|)) would yield a negative entropy 
at the Unruh- Hawking temperature, which is very hard to accept on a physical ground. 
Summarizing, it seems to us that this second approach, which is the natural generalization 
of the procedure used in the scalar case, does not yield a correct result. 



IV. SUMMARY AND DISCUSSION 



The main result of this paper is the proof that the optical method (the "first approach" ) 
can be used to compute one-loop quantities in the Rindler space also in the case of the 
photon field. The method has been developed employing a general covariant gauge choice. 
Furthermore, by a comparison with other methods, we have seen that this method produces 
the same result as the point-splitting procedure. 
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It is also important to stress that the partition function arising from our method is com- 
pletely free from "Kabat's" surface terms. This is very important because, as we previously 
said, the approaches based on the direct computation in the Euclidean Rindler space us- 
ing ^-function or heat-kernel techniques produces such anomalous terms and further 
regularization procedures seem to be necessary to get physically acceptable results. 

We have also developed a general optical formalism for the Maxwell field in the covariant 
gauges based on Hodge de Rham formalism which, in principle, can be used in different 
manifolds than the Rindler space. 

However, many problems remain to be explained. In particular, both in the photon 
and in the scalar case the relation between the optical approach and the direct approach in 
the manifold with the conical singularity remains quite obscure. This is due to difficulties 
involved in computing the Jacobian of the conformal transformation in the presence of 
conical singularities. Moreover, while the optical approach can be used in the case of massless 
fields without particular difficulties, as soon as the fields have a mass the optical method 
becomes much harder to apply. In this case, the direct computation in the manifold with 
conical singularities could show its advantages, provided one knows how to compute the 
above Jacobian. 

Another general point which requires further investigation is the request of self- 
consistency of the thermodynamics of the gas of Rindler particle, when the temperature 
is not the Unruh one. This is a very important point in calculating the correction to the 
entropy of a black hole supposing such corrections due to the fields propagating around 
it. We remind one that the Rindler metric approximates the region near the horizon of a 
Schwarzschild black hole. The entropy of the fields is computed using the relation (where (3h 
is the Unruh-Hawking temperature, 2n in the Rindler case): Sp H = f3\dpFp\p H . In calcu- 
lating the previous derivative at (3 — j3jj, one has to consider also the partition function off 
shell, namely evaluated at (3 7^ 2tt and f3 near (3h- It is not so clear whether it is necessary 
or not that the thermodynamical laws hold also for (3 7^ (3h and (3 near (3h in order to assure 
the consistency the procedure followed in calculating the entropy of the fields at (3 = /3#. 
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Moreover, it is well known that the off shell quantum states of a field are affected by several 
pathologies on the horizon event .[] Furthermore, they are unstable states in a semiclassical 
approach to quantum gravity due to the divergence of the renormalized stress tensor on the 
horizon. Thus, it is reasonable to wonder about the thermodynamical consistency of the 
results when one works off shell. We conclude with a discussion on this point. 

Let us first consider the point-splitting result and the annoying dependence on the pa- 
rameter £ of the massless scalar field results and its relation with the request of a consistent 
thermodynamics. As we have already said, on more regular manifolds the (integrated) phys- 
ical quantities should not depend on the actual value of £, whereas in the case of Euclidean 
Rindler wedge the conical singularity introduces an, apparently unphysical, £ dependence 
in the physical quantities.^ A similar problem occurs even in flat spaces in the presence 



of boundaries ||21|| . In those cases, one can see that the renormalized energy-momentum 
tensor diverges on the boundary, unless the coupling is conformal. So, one could say that 
the conformal coupling is, in this sense, more "physical" than the others. By this we mean 
that it behaves like a real field, such as the Maxwell field. 

Inspired by this fact, we shall look for a criterion to choose a value of £ which is more 
"physical" than the others as far as the thermodynamics is concerned. Thus we shall discuss 
the consistency of the thermodynamics of the point-splitting results. 

From the thermodynamics, we know that we can obtain the internal energy Up'~g' (see 
Eq. (fj)) as the derivative with respect to (3 of an appropriate free energy F^ s ' multiplied 



7 Rindler thermal states with (3^2^ violates several axioms of the QFT in curved backgrounds. 



For example, see 1 33 1 and ref.s therein. 

8 Notice that the case of the cosmic string theory is quite different because different values of £ 
correspond to different internal structures of the string. This is obvious by considering a string 
with a finite thickness, which has a nonvanishing curvature within itself. In the limit of a vanishing 
thickness, the curvature R gets a Dirac 5 behaviour along the string in the Lorentzian manifold. 
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by ft, possibly corrected by an suitable energy-subtraction procedure. Taking the space 
homogeneity along the y and z directions into account, we found the form of the free energy 



pp.-s 



LyL Z 



28807r 2 e 2 

LyL Z 

28807r 2 e 2 



j) 



'2V 

J. 



30(6£- 1) + 3 



W(£,e) + 



/3 



(50) 



The unknown function /(£, e) can be dropped by requiring that the entropy Sp^ = ft 2 dpF^ s ' 
vanishes at ft — > +oo. The function U, which does not depend on ft but can depend on 
the geometry background, is necessary due to the fact that the energy in Eq. is the 
Minkowski renormalized one, but we want to remain on a more general ground in order to 
use the thermodynamical laws. In other words, we may notice that the energy in Eq. 
becomes negative if the temperature is sufficiently low, for example, in the most interesting 
range < £ < 1/6, and hence such an energy cannot directly arise form a statistical partition 
function but a further subtraction procedure must have taken place. The function IA takes 
into account this energy subtraction procedure. 

From statistical thermodynamical laws, one expects that the y and z principal pressure, 
namely T yy and T zz in Eq. ([[]), integrated over dzdr^Jg and dydr^fg respectively, can be 
obtained taking the L y (L z ) derivative of the previous free energy, with the sign changed 
and U opportunely chosen. An easy computation shows that, due to the terms containing 
ft~ 2 , this does not hold for any value of £, but only in the conformally coupled case, £ = 1/6. 
After the Minkowskian energy subtraction, the corresponding free energy reads 

A i 



rpp.s. 



^.-s. sub 



'27rV 

J +3 



(51) 
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This is just the free energy obtained by the optical method after the Minkowski renormaliza- 
tion. Therefore, it seem that only in the conformally coupled case the stress tensor (|1]) yields 
a consistent thermodynamics, at least as far as the relation between energy and pressures is 
concerned. 

Now, let us consider the photon case. In such a case we have not the freedom to adjust a 
parameter in the stress tensor in order to agree with the thermodynamics. The free energy 
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we find from the total energy in the case of the photon stress tensor of Eq. flT5|) reads: 



^hot. p.-s. 

P 
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U{e) + M 



(52) 



As before, we can drop the term containing the undetermined function /(e) by requiring a 
vanishing entropy in the limit (5 — * +oo. The above free energy produces the point-splitting 
internal energy and, after the Minkowski renormalization, it coincides with the free energy 
obtained by renormalizing that obtained by the optical approach, Eq. (|43|) . 

The point is that if we apply the above procedure to compute the integrated principal 
pressures along the y and z directions to the above photon free energy, there is no way to 
choose U in such a way to get the same result as integrating the yy and zz components of 
the photon stress tensor in Eq. fllS]). This is due to the presence of a term proportional to 
f3~ 2 and the independence on (3 of the function U. 

In order to get the "correct" pressures (but a wrong internal energy!) employing the 
derivatives as previously pointed out, one should take a free energy which is twice that in 
Eq. fl5"0|) with £ = 1/9, f = and U opportunely chosen. 

Hence, it seems that the point-splitting stress tensor of photons in the Rindler wedge 
does not give a consistent thermodynamics.^ It is very important to remark that the above 
thermodynamical argument cannot be applied to the cosmic string theory, since in that case 
the stress tensor in Eq. ( [Tlf) is the zero-temperature one, and j3 is not the inverse of the 
temperature. 

In a pessimistic view, this problem and the £ dependence of the integrated quantities 
in the scalar case could be considered as another proof of the inconsistency of the Rindler 

9 This problem arises also dealing with the massless spinorial field as it simply follows from the 
point-splitting renormalized stress tensor obtained in || (analytically continued from the cosmic 
string to the Rindler space). 
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theory (and maybe of the Schwarzschild theory) when one works at temperatures different 
from the Unruh-Hawking one, and a discouraging result for the attempt to evaluate the 
correction to the Bekenstein-Hawking entropy through the "off shell" procedure. 10 
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V. APPENDIX 

In computing the photon ( function on S 1 x if 3 one meets the ( function of a scalar field 
in the same background, both in conformal and minimal coupling. Therefore, it is useful 
to report here these ( functions. The small fluctuations operator for a scalar field in the 
optical metric is 

Lt = A - 6£ = -[d 2 T - rd r + r 2 d 2 r + 6f]. 

where A is the Hodge de Rham Laplacian on S 1 x H 3 . A complete set of eigenf unctions has 
been given in the main text, Eq. (|35|) , with eigenvalue \u\ + u 2 + 1 — 6£]. Therefore, the 
local ( function is 



00 roo r 

((s\Le)(x)= ]T / duo j d 2 k[u 2 + uj 2 + l-6^- s 4>*(x)4>(x) 

n=-oo ^° J 

r(s) ^.^JO 



10 However, it could be possible to interpret the entropy formula at Hawking's temperature, with- 
out making use of thermodynamical laws off shell. Maybe, possible ways could arise studying the 
geometrical entropy employing the replica trick i"0|,22|. 
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7T r(s-|) /27r\ 3 - 2S 



3 /3 r — \ / /3 " 2 "' 



8tt 2 /? r(s) V £ / 

where E^s; a) = I]^L [ n2 + fl2 ]~ s is the Epstein £ function. In the conformally coupled case, 
the Epstein function becomes a Riemann ( function and so 



C-fo*) s CW^)<*> = ^(!) " ^0.(2.-3), 
One can easily check that ( c - c - s -(s; x) | s=0 = and 

d 7T 2 



Another important case is the minimally coupled one, £ = 0, for which there is not a more 
explicit form. However, using the identity 

and the fact that the MacDonald function K v {x) is analytic in the index v and decays 
exponentially as \x\ — > oo so that the third term in the previous expansion is analytic in s 
(and vanishes as s — > 0), we find that the £ function does not vanish in s = 0: 

We do not know the value in zero of the derivative, but it is not required in our computations. 
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